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Based on the Kotliar-Ruckenstein slave-boson scheme we develop a configuration-interaction (CI) 
approach which is suitable to improve the energy of symmetry-broken saddle-point solutions. The 
theory is applied to spin-polaron states in the Hubbard model and compared with analogous results 
obtained within the Hartree-Fock approximation. In addition we show that within the infinite D 
prescription of the Gutzwiller method a CI approach does not improve the variational result since 
in the thermodynamic limit matrix elements between different inhomogeneous states vanish due to 
an 'orthogonality catastrophe'. 
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I. INTRODUCTION 

The Gutzwiller Ansatz is a variational wave function 
for correlated electronic models with purely local inter- 
actiontiiS The basic idea to treat these Hubbard-type 
hamiltonians is to partially project out configurations 
with doubly-occupied sites from the Fermi sea in order 
to optimize the contributions from kinetic and poten- 
tial energy. As a consequence, in contrast to the con- 
ventional Hartree-Fock (HF) theory, the Gutzwiller wave 
function captures correlation effects like the band nar- 
rowing already on the variational level. However, the 
exact evaluation of the ground state energy within the 
Gutzwiller wave function is fairly difficult and up to now 
has only been achieved in one and infinite dimensions.'^ In 
the latter case the solution is equivalent to the so-called 
Gutzwiller approximation (GA) which has been applied 
to describe a variety of finite dimensional systems rang- 
ing from the properties of normal '^He (cf. Ref. [1) to the 
stripe phase of high-Tc cuprates^iS, 

The GA in its original formulation was restricted to ho- 
mogeneous paramagnetic systems and only later on gen- 
eralized to arbitrary Slater determinants by Gebhard? 
and, more recently, by Attaccalite and Fabrizio.? The 
same energy functional was obtained from the Kotliar- 
Ruckenstein (KR) slave-boson formulation of the Hub- 
bard model when the bosons are replaced by their mean- 
values.— Unconstrained minimization of the KR (or 
Gebhards) energy functional on finite clusters in gen- 
eral yields inhomogeneous solutions which break transla- 
tional and spin-rotational invariance.^'^*^ This approach 
has been used for the investigation of electronic inhomo- 
geneities, such as stripes and checkerboardsi^ii^iii, in the 
context of high-Tc superconductors. 

Incorporation of fiuctuations in the frame of the time- 
dependent Gutzwiller approximation tends to restore the 
original symmetry of the systemJ^ An alternative would 
be the construction of a wave-function which is a lin- 
ear superposition of equivalent symmetry- broken states. 
In case of stripe states^^'^^ one could e.g. envisage a 
superposition of solutions which are translated perpen- 
dicular to the stripe direction and also the correspond- 
ing solutions which are rotated by 90 degrees. In case 



of the unrestricted Hartree-Fock approximation such a 
configuration-interaction (CI) method has been proposed 
in Ref. [l^ and applied to the case of stripe textures in 
Ref. [O. 

The present paper investigates the possibility wether 
an improvement of the inhomogeneous Gutzwiller ap- 
proximation is possible within an analogous framework. 
In Sec. II we evaluate the matrix elements of the Hub- 
bard hamiltonian between different inhomogeneous solu- 
tions obtained from the saddle-point approximation of 
the KR slave-boson scheme. Based on these results we 
construct a a CI method which in Sec. HI is applied 
to spin polaron states. We compare ground state ener- 
gies with exact diagonalization results and for larger lat- 
tices evaluate the dispersion relation of the spin polaron 
states which can be compared with analogous solutions 
obtained in the tJ-model. In this context we also compare 
our results with angle-resolved photoemission (ARPES) 
experiments on Sr2Cu02Cl2. 

In appendix [Xl it is shown that the infinite D prescrip- 
tion of the Gutzwiller approximatioiti cannot be used for 
an analogous construction of a CI approach. The reason 
is that in the thermodynamic limit this scheme leads to 
an 'orthogonality catastrophe'— so that energy correc- 
tions and the dispersion of quasiparticles vanish. 



II. MODEL AND FORMALISM 

Our investigations are based on the one-band Hubbard 
model 

ij,a i 

where c^^l destroys (creates) an electron with spin cr at 

site i, and rii o- = cJ^Ci.o-. U is the on-site Hubbard 
repulsion. 

Following KRii. we enlarge the original Hilbert space 
by introducing four subsidiary boson fields e[^\ s^^^ , s^-^^', 

and d'l'^ for each site i. These operators stand for the 
annihilation (creation) of empty, singly occupied states 
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with spin up or down, and doubly occupied sites, respec- 
tively. Since there are only four possible states per site, 
these boson projection operators must satisfy the com- 
pleteness condition 



Furthermore 



4^s^<T + did. 



(2) 



(3) 



Then, in the physical subspace defined by Eqs. (|2l3p the 
Hamiltonian ([T]) takes the form 



H 



(4) 



The expectation values of the z-operators Eq. 



(10) 



are equivalent to the renormalization factors derived 
within the infinite V prescription of the Gutzwiller 
approximation^ (cf. Eq. (jA10[) in appendix E)) . 

In previous works^iiS we have proposed a method for 
minimizing the KR energy functional E" = 
on finite clusters without imposing constraints with re- 
spect to translational and spin rotational invariance. In 
the remainder of this section we evaluate the matrix el- 
ements of the Hubbard model between two different in- 
homogeneous solutions 1^1/") which then will be used in 
order to partially restore these symmetries. 

We start with the overlap between wave- functions be- 
longing to different inhomogeneous solutions 



with 



(5) 



and has the same matrix elements than those calculated 
for ^ in the original Hilbert space. The operators f^^^ 
are the electron annihilation (creation) operators in the 
new Hilbert space. 

In the saddle-point approximation we can represent the 
wave-function for a specific inhomogeneous solution a as 



(6) 



where |$q) is a Slater determinant and the bosonic part 
\Bq) is a coherent state 

\B^) = eE,('Jr4+E.^t.<„+s?^^-i/2)|o). (7) 



= ($o"|<)(So"|5o^) 
where the overlap between coherent states reads as 



(11) 
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The fcrmionic overlap is given by 

($o"|$^) = (<i>ol*^>T('i>ol*[^> 



(12) 



(13) 



and the evaluation of the spin-dependent factors is out- 
lined in appendix [Bl 

We now proceed by calculating the matrix elements of 
the hamiltonian Eq. ^ in the basis of the inhomoge- 
neous wave- functions l^*"). From the above definitions 
one obtains for the Hubbard interaction 



Since a coherent state contains an arbitrary number of 
bosons the constraints Eq. (|2l3p are only fulfilled on av- 
erage for a given inhomogeneous solution a provided that 

1 = {e?r + J2(~'Z'^)' + (dff 

(T 

Note that here and in the following expectation values of 
fermion operators are denoted with respect to the Slater 
determinant of /-electron operators. 

The problem with the Ansatz Eq. ([6|) is that one does 
not recover the correct non-interacting limit [/ — > for 
which — > 1. Therefore KR^^ introduced a unitary 
transformation in order to represent the z-operators in 
Eq. dll) as 



1 



= (e|si,<,-t-4-„rfj) 



1 



did, + sl„s.,^a 



(8) 



so that 



(14) 



The kinetic term is evaluated in a similar way as 



with the fermionic part 



('i>olc. 



'^l.cr'^j,^ 



(15) 



(16) 



and the brackets are defined in Eq. jBTj in appendix [B] 
The matrix elements of the 'bare' bosonic 'z'-operators 
from Eq. ([5]) read as 



q/3 



dfs: 



(17) 



Now we have to deal again with the problem that the 
z-factors as defined in Eqs. P?|) do not yield the uncorre- 
lated limit, i.e. = zff — > 1 for [7^0. It is straight- 
forward to proof that the representation of Eq. ([5]) does 
not work in this case since the above limit is only obeyed 
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for homogeneous paramagnetic solutions. However, due 
to a non-symmetric population of momentum states on 
finite clusters or in case of inclusion of an electron-phonon 
coupling the charge and spin structure in general is in- 



J 



homogeneous even in the limit f/ — > 0. 

A possible representation which yields 
U is given by 



1 for 



1 ^1*^4 ^i,-a^'i,-a 



1 - 44 - 4,a4.a 4-''h-'r 



1 



I 



(18) 



Note that in the physical subspace defined by Eq. ([2]) 
the square root factors are identically 'one'. On the other 
hand, upon evaluating the matrix elements of Eq. (jlSp 
between coherent states a, (3 one obtains the hopping 
renormalization factors 



a/3 



{BS\4,M 



1 



(19) 



(df)2((n) 



1 - {rit-a 



1 - (ni 



(1-Mf + (rfT)')((")".-(df)2) 



where we have used the constraints Eqs. (|2I3|) to replace 
the boson fields but df by fermionic expectation values. 
The 'z-factors' Eq. (I19|l show the correct behavior z"^ 
1 for [/ — > and the diagonal elements reduce to the KR 
renormalization factors Eq. PH)) . i.e. zf" = gf^.. 

In appendix[X]it is shown that the renormalization fac- 
tors Eq. pop can be also motivated from the generalized 
Gutzwiller approach in the limit 2? — > oo. 



III. RESULTS 

In the previous section we have calculated the matrix 
elements between different inhomogeneous states j^*") of 
the Hubbard model. These results are now used for evalu- 
ating an improved ground state energy and wave- function 
similar than in the configuration interaction approach 
based on unrestricted HE wave- functions. 

We apply the method to the investigation of spin po- 
laron states on a square lattice, i.e. we have one hole with 
respect to half-filling. Minimization of the KR (or GA) 
energy functional leads to the localization of this hole at 
a given site Ra (cf. Ref . |^ for a method of performing the 
unrestricted variation) and we denote the corresponding 
projected or fermion-boson wave-function with l^*"). 

Now we generate all translations of this solution within 
the same sublattice since solutions belonging to different 



sublattices are orthogonal. The superposition 



(20) 



thus only includes states j^*") with the same energy 
E = E°' . In principle one could systematically improve 
the approach by including also excited states of the un- 
derlying fermionic Slater determinant. 

If we apply the hamiltonian Eq. ((T|) to Eq. (^0)) one 
obtains the following eigenvalue problem 

{■^^\H\^^)vp=eS^pvp (21) 

where the matrix 5*^^ is defined in Eq. (|lip . 

A. One hole states in the 4x4 lattice 

We start by investigating the quality of the present 
approach with regard to exact results and the HF con- 
figuration interaction method (CIHF). 

Table IIII Al reports the energy correction obtained 
with our slave-boson configuration interaction approach 
(CISB) as compared to the unrestricted GA. The values 
for the exact result, the CIHF and the unrestricted HF 
(from Ref. TJ) are also shown for comparison. 



U/t 


exact 


HF 


GA 


CIHF 


CISB 


4 


-0.91658 


-0.83139 


-0.88815 


-0.83501 


-0.89091 


6 


-0.74794 


-0.64222 


-0.70020 


-0.66214 


-0.70497 


8 


-0.634203 


-0.52884 


-0.57518 


-0.54767 


-0.60295 


16 


-0.42546 


-0.33589 


-0.37130 


-0.34604 


-0.38091 


32 


-0.308473 


-0.23160 


-0.27209 


-0.23627 


-0.27685 


50 


-0.266039 


-0.19335 


-0.23954 


-0.19617 


-0.24362 



TABLE I: Energy per site for 15 particles on a 4 x 4 lattice. 
The values of the exact result, HF and CIHF method have 
been taken from Ref. [l^ 

It turns out that the CISB leads to an energy cor- 
rection to the GA result which is of the same order of 
magnitude than the CIHF correction to the HF energy. 



4 



O O CIHF (U/t=10) 
° •-*CISB (U/t=10) 

Q-B CISB (U/t=20) P 
O A CISB (U/t=40) 



.0) 



« o o o e o •' 

A-,. ..-A 
* A^A -A A'Ar * 



(7C,0) 



A. A^_ 



B' - 



(7t,7t) (lt/2,7C/2) (0,0) 



FIG. 1: (color online). Dispersion of the spin polaron in the 
Hubbard model evaluated within the CISB {U/t = 10, 20, 40) 
and CIHF {U /t = 10) method. Energies are with respect to 
the half- filled antiferromagnet. 



However, this improvement is on top of the GA which it- 
self provides a much better estimate for the ground state 
energy than the HF approximation. For example, one 
finds that for U/t — 8 the CISB differs from the exact 
result by ~ 5% whereas it is « 13% in case of the CIHF. 



B. One hole states in the 16x16 lattice 

We continue by evaluating the dispersion of the 
spin polaron in a 16 x 16 lattice. This prob- 
lem has been extensively investigated within the tJ 
model,^^^^^>25^^^ where for small J/t one finds 
a bandwidth ~ J which turns over into a 2t^ / behavior 
for large J/t. Furtheron the dispersion is characterized 
by a maximum at (0, 0) (and the analogous (tt, tt) point) 
and displays a 'hole pocket' at (7r/2, 7r/2) which is slightly 
lower in energy than the (tt, 0) point. 

Fig. [T] displays the polaron dispersion obtained within 
the SBCI method for U/t = 10,20,40. For compari- 
son we also show the U/t = 10 result obtained from 
the CIHF method. Since the wave-function incorpo- 
rates only polaron states localized on the same sublattice 
the dominant contribution to the dispersion is given by 
Ek ~ 'it' cos{kx)cos{ky) + 2t[cos{2kx) + cos{2ky)]. There- 
fore at the point k = (7r,7r/2) the energy difference be- 
tween CIHF and CISB corresponds to the difference be- 
tween GA and HF energies for the spin polaron. Since 
within the CISB approach the matrix elements which en- 
ter Eq. (|2ip are additionally scaled by the bosonic ex- 
ponential overlap Eq. (fT2|) the corresponding long range 
contributions to the dispersion are in generally smaller 
than for the CIHF method. On the other hand, this scal- 
ing affects also the matrix Sa0 in Eq. pTjl so that due to 
partial cancellation the overall effect on the bandwidth 
is less pronounced as one might expect (see below). 

From analogous investigations in the tJ- 
model2a^>22^>2i^>26^^ it is known that the 



dispersion of a single hole has a saddle-point at 
k = (tTjO) and k = (7r/2,7r/2), where the latter corre- 
sponds to the minimum of the band. From Fig. [T] it 
turns out that the CIHF spin-polaron dispersion also 
displays the minimum at fc = (±7r/2, ±7r/2) whereas 
within the CISB method the state at fc = (±7r, 0), (0, ±7r) 
is slightly lower in energy. However, a direct comparison 
of results between tJ- and Hubbard model is hampered 
by the fact that the strong coupling expansion of the 
Hubbard model generates a three-site term of order J in 
addition to the 'conventional' tJ-model. Since we find 
that the energy difference between fc = (±7r/2, ±7r/2) 
and fc — (±7r, 0), (0, ±7r) states is always smaller than 
J = At'^/U there appears no inconsistency with results 
from the tJ-model. In fact, calculations of a single hole 
in the antiferromagnet based on an expanded tJ-model 
(including the three-site term) provide evidence that the 
minimum of the band may be at fc = (±7r,0), (0,±7r)— . 
This finding is also substantiated by exact diagonal- 
ization results of the same model on small clusters ^ 
Unfortunately, for the full Hubbard model there are no 
conclusive answers from Quantum Monte Carlo or exact 
methods yet availablcj ^^'^^ 





SCBA 


CIHF 


CISB 


J 


E(tt/2,tv/2) 


W 


-E(^/2,7r/2) 


W 


E(tv/2,-!t/'2) 


W 


0.1 


-2.785 


0.239 


-1.84 


0.231 


-2.4786 


0.263 


0.2 


-2.540 


0.430 


-1.703 


0.513 


-2.204 


0.421 


0.3 


-2.360 


0.600 


-1.588 


0.817 


-2.036 


0.68 


0.4 


-2.209 


0.741 


-1.487 


1.118 


-1.95 


1.031 



TABLE II: Binding energy Epoiaron — Eaf taken at momen- 
tum q = (7r/2,7r/2) and the bandwidth W for various values 
of J = 4t'^/U. Shown are results for the self-consistent Born 
approximation (SCBA) of the tj-model (from Ref. [20l) and 
the CIHF and CISB method for the Hubbard model, respec- 
tively. 



Table [III reports the bandwidth, and the energy at 
fc = (±7r/2, ±7r/2) of the spin polaron dispersion ob- 
tained within the SCBA,Sfl CIHF and CISB method, re- 
spectively. Note that for the latter approach the band- 



width isW = E, 



(0,0) 



E, 



(tt.O) 



whereas for the SCBA and 



E(tv/2,q)- De- 



CIHF methods it is given by W = i?(o,o) 
spite this difference we find that the CISB bandwidth 
scales as ly « 2.2 J up to J « 0.3 in agreement with anal- 
ogous considerations in the tJ-model. It also turns out 
that (at least for J > 0.1) the CISB bandwidth is smaller 
than that of the CIHF approach. Formally this is again 
due to the additional renormalization of the matrix ele- 
ments by the bosonic exponential overlap Eq. p2| . On 
the other hand it is quite natural that the CISB approach 
leads to 'heavier' spin polarons than the CIHF method 
due to the incorporation of correlation effects already on 
the Gutzwiller level. Similar to the case of the 4x4 
lattice the CISB leads to a significant energy correction 
with regard to the CIHF as exemplified by the value of 
E{tt/2,7v/2) in table [ni 
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C. Comparison with experiment 

Undoped cuprate superconductors are antiferromag- 
netic Mott insulators. Within a angle-resolved photoe- 
mission (ARPES) experiment, one can in principle ob- 
serve the dispersion of the created hole in the antiferro- 
magnetic background of these compounds and compare 
with that of of the spin polaron quasiparticle concept 
from the previous section. On the basis of the single- 
band description it is now well established from LDA^^ 
and the analysis of ARPES data^** that a next-nearest 
neighbor hopping t' has to be considered in the model. 
In particular, it has been found^** that the quasiparticle 
dispersion from (7r,0) to (7r/2,7r/2), which is determined 
by t' , is characteristic for the different cuprate families. 
Our analysis below is therefore based on the extended 
Hubbard model, which corresponds to Eq. ([T]) when the 
hopping tij is restricted to nearest t and next-nearest 
~ t' hopping. In Fig. [5] we fit the resulting spin polaron 
dispersion to ARPES data on undoped Sr2Cu02Cl2 ob- 
tained Wells et ai. Since the experiment measures the 
single particle Green's function for electrons the disper- 
sion in Fig. [5] is 'reversed' with respect to those shown 
in Fig. [1] which were obtained for holes. 

We can use the experimental energy differences AEi — 
E(n/2,7T/2)~E(ofi) and AE2 = -B(7r/2,7r/2) -^^(7r,o) in order 
to fit two of the three parameters (t, t' , U) . Therefore we 
additionally use our results from Ref. |36| where we have 
fitted the magnon dispersion of undoped La2Cu04 within 
the time-dependent Gutzwiller approximation. In this 
case the value of the Hubbard repulsion U/t ^ 8 could 
be accurately determined from the dispersion of spin ex- 
citations along the magnetic Brillouin zone whereas this 
dispersion is rather unsensitive to t' . Given that the Cu 
onsite repulsion should not depend very much on the ma- 
terial we there also use the ratio U/t in our present fit 
of the spin polaron dispersion for Sr2Cu02Cl2. As a re- 
sult we find that the ratio t' /t = —0.2 yields an overall 
good agreement with the data and the nearest neighbor 
hopping t = SOOmeV is set by the absolute energy scale. 
The ARPES data in addition allow for an accurate deter- 
mination of t' so that a combination of both approaches 
in principle can be used to obtain parameter sets for the 
Hubbard model in order to describe different materials. 



IV. CONCLUSIONS 

We have developed a configuration interaction ap- 
proach based on the KR slave-boson mean-field formu- 
lation of the Hubbard model^^. In principle this method 
provides a controlled scheme for including fluctuations 
beyond the mean-field solution. Formally this has been 
achieved by several authors within the functional integral 
formalism.— i^Si^i^i^ Here we have discussed an alter- 
native extension which is based on the observation that 
unrestricted variation of the KR energy functional in gen- 
eral leads to a class of degenerate solutions which are con- 
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FIG. 2: (color online). Dispersion of the spin polaron in the 
extended Hubbard model evaluated within the CISB {U/t = 
8, t'/t = -0.2, t = SOOmeV). The right panel shows the 
direction along the boundary of the magnetic Brillouin zone. 
Experimental data are from Ref. [ssl . 



nected by symmetry transformations. The CISB method 
discussed in this paper allows for a tunneling between 
these degenerate solutions and thus for a construction of 
eigenstates with well defined momentum. 

Although the KR mean-field energy functional is iden- 
tical to the that obtained with the generalized Gutzwiller 
wave- function in P — > oc^ the considerations in appendix 
lAl show that the latter approach leads to an 'orthogonal- 
ity catastrophe' for matrix elements between different in- 
homogeneous states. Therefore one would have to invoke 
1 I'D corrections in order to construct a CI approach also 
within the Gutzwiller method. 

Application of the CISB to the spin-polaron problem 
for the Hubbard model leads to a significant energy gain 
with respect to the CIHF method. In addition we have 
obtained a minimum of the spin polaron dispersion at 
k = (±7r, 0), (0, ±7r) in contrast to analogous calculations 
in the t J-model but also in contrast to the CIHF method. 
However, calculations based on the full strong coupling 
expansion of the Hubbard model ,^^1'^° which take into 
account the three-site terms of order t^ /U, neglected in 
the 'conventional' tJ-model, indicate the occurence of 
dispersion minima around the corners of the magnetic 
Brillouin zone. To our knowledge there are no recent ex- 
act diagonalization studies of one hole in a \/l8 x \/T8 
or X \/20 Hubbard cluster which could substantiate 
the findings of Ref. [s^- However, since on the mean- 
field level the KR slave-boson formulation of the Hub- 
bard model takes into account correlations beyond HF 
we expect that the CISB is more accurate concerning fine 
details of the spin polaron dispersion as compared to the 
CIHF method. Further investigations are needed in or- 
der to confirm the finding of one hole dispersion minima 
at fc = (±7r,0), (0,±7r) in the Hubbard model. 

Finally, we have included a next-nearest neighbor hop- 
ping t'/t < in the bare hamiltonian in order to fit the 
low enery dispersion of Sr2Cu02Cl2 from ARPES exper- 
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imentsj^ The parameter t' is essential in order to obtain 
the measured dispersion along the border of the magnetic 
Brillouin zone. More recent ARPES experiments^ have 
also revealed a strong dispersion along the (0, 0) — + (tt, 0) 
direction. Within a one-band description modeling of 
these data requires inclusion of a significant third near- 
est neighbor hopping. However, since our CISB approach 
can be implemented also on the more realistic three-band 
model it would be interesting to study the spin polaron 
dispersion within this hamiltonian. The comparison with 
ARPES experiments would then allow to elucidate the 
parameters of this hamiltonian for different cuprate ma- 
terials. Moreover, since the superposition in Eq. ([20]) 
can be extended to include also excited states, it should 
be possible to calculate also the incoherent part of the 
ARPES spectrum and thus to provide a more detailed 
description of the data. Work in this direction is in 
progress. 
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APPENDIX A: GENERALIZED GUTZWILLER 
APPROXIMATION 

Following Ref . the Ansatz for a given inhomogcncous 
state a can be written as 



BT 



9 



9 



(Al) 



(A2) 



where the uncorrelated state |<i>o ) is a Slater-determinant 
with an inhomogeneous density matrix a and Di — 
rii^ini^l is the double occupancy operator. For later pur- 
poses we also define the operators for single occupied 
(with spin a) and empty sites: 



(A3) 
(A4) 



The parameters ^f^ and r]f have to be determined vari- 
ationally. Gebhard^ has shown that the requirement 
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2K{a) ^ 



\ln [n-a;f(A 



(A5) 



leads to the same energy functional than the Kotliar- 
Ruckenstein slave-boson approach in the mean-field ap- 
proximation when the expectation values are formally 
evaluated in the limit of infinite dimensions. Here D^^'" 
denotes the Hartree-Fock decoupled double occupancy 
operator in the basis of the Slater determinant \^q). Eq. 



[5]) yields a relation between the variational parame- 
ters (7, nf„, rjf and the variables xf which turn out to be 
the relevant parameters when one evaluates expectation 
values in infinite dimensions. The essential step in this 
direction is to express the operator Bf defined in Eq. 
i2\\ in terms of the xf as 



(A6) 



E,A l + x?(D, 



and the expectation values are defined with regard to 
|$q). An important result of the d — *■ oo description 
is the equivalence of local densities in the projected and 
unprojected states 



which will be used in the following. 

First the double occupancy can be evaluated as 

(vI/"|A|*"> = -Df - (A)" (1 + xf{E^r) 



(A7) 



(A8) 



which allows one to perform the variations with respect 
to the double occupancy Vf instead oixf (or g, ^if^, and 
ryf ). Analogously the hopping term of Eq. ([l} is given 
by 

m4,.<^3..m - <it.<iiAn\4,.c,An) (A9) 

with the hopping renormalization factors 



1 



(AlO) 



Similar than in Eq. (jA8[) expectation values of a pro- 
jection operator Pi = Di, Si,a, Ei with regard to |^") 
have been denoted with calligraphic letters. The above 
representation of the hopping factors allows for a inter- 
pretation of the renormalized kinetic energy in terms 
of 'probability ratios'4ii2, Consider the term qf„Cicr\^Q) 
which is the sum of two processes: The contribution 

Cia{l — (rii^-a)") originates from the annihilation of 
a singly occupied (and thus creation of an empty) site 
and is weighted by the ratios between projected and un- 
projected probabilities of this process. The contribution 
^ Cio-(rii__cr)" weights in a similar way the annihilation 
of an electron on a doubly occupied site. 

We now proceed by evaluating the matrix S which con- 
tains the overlap elements of wave- functions belonging to 
different inhomogeneous states 

Sa.p - (vl/^lM^^) =n(<&o|Sfsfl*o) (All) 
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+ 



.13 



(<folA|< 



where we have used Eqs. (|A6|A7IA8P and the fact that 
only local contractions survive in infinite dimensions. Eq. 
(ElT|) also requires the evaluation of matrix elements of 



Pi between different Slater determinants (^q |Pi|$o) 



For 



example, one finds for the double occupancy operator 

($^|A|<f(?) = [r^.,T]„^[na]«^ (A12) 



and the brackets are defined in Eq. (|B7p . 

Schwartz's inequality together with the relation be- 
tween harmonic and geometric mean 



yields 



(A13) 
(A14) 

(A15) 



where the equals sign holds for a — (3. 

Analogously to S one can evaluate the matrix elements 
of the Hubbard hamiltonian Eq. ([T]). For the double 
occupancy operator one obtains 



0/ r. 
- — ^aP 



(A16) 



and the matrix elements of the hopping term are given 

by 



In principle it is possible to evaluate the matrix el- 
ements from Eqs. (|A17p in terms of the Slater de- 
terminants |$o)j however, the calculation of contribu- 
tions which involve density correlations of the form 
($Q |ni,_o-nj,_o-|$o) are rather time consuming. We 
therefore simplify the expression of the projections Eqs. 
(|A18IA19[) by the following argument. With regard 
to the matrix element Eq. (jA17[) the projection Eq. 
(lAlSp describes the annihilation of a particle with spin 
a in the Slater determinant ($ol- The two contribu- 
tions measure the probabilty wether site i in the state 
a is singly or doubly occupied. Accordingly we replace 
the corresponding projections by their mean- values, e.g. 
1 — rii^^a 1 — {rii-a)"'- In the same way Eq. (|A19|) 
describes the annihilation of a particle with spin a in 
the Slater-determinant I'I'o)'' and we approximate in this 
case l — rii^^u ^1— {ni^^„)^ . Within this approximation 
one obtains for the projected creation and annihilation 
operators 



a/3 t 
Ba 



(A20) 
(A21) 



where the q"^ are equivalent to the renormalization fac- 
tors Eq. (fTO)) derived with the KR slave-boson method. 

In case of the GA we observe from Eq. (jAlip that 
Sa^p is a product over lattice sites of terms less than 
'one' which in the thermodynamic limit leads to an 'or- 
thogonality catastrophe'^^ and thus Sap = Sap- There- 
fore we find that within the 'infinite D' prescription of 
the Gutzwiller approximation^ different inhomogeneous 
states are orthogonal to each other. As a consequence it 
turns out from Eqs. (jA16IA17|) that these states are not 
connected by matrix elements of the Hubbard hamilto- 
nian so that a CI approach does not yields any correction 
to the symmetry-broken solutions. 



(A17) 

Using Eqs. (|A6IA7IA8p the projections of the creation 
and annihilation operators can be expressed as 



B'^c'. B^ 



Bfc,„B 



/3 



(1 



CQ CP 



(5,,)"(i?,)/3 



(1 



CO. QP 



(A18) 



(A19) 



APPENDIX B: FERMIONIC MATRIX 
ELEMENTS 

When we restrict to collinear inhomogeneous 
Gutzwiller solutions, i.e. where the associated den- 
sity matrix is diagonal in spin space, we can represent 
the non- interacting state |<I>q) as 



Q.t CK.t Ct.t 

ai;>2,>3,.T • 



''Na. 



JO) 



(Bl) 
(B2) 



and the operators ^ are related to the real space op- 
erators Ci^a- by the linear transformation 



E 



(B3) 



which defines the specific inhomogeneous solution. De- 
tails for the calculation of the amplitudes 4>f^{k) within 
the Gutzwiller approximation can be found in Ref. [^. 
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Within these definitions the evaluation of matrix ele- 
ments between different Slater determinants is analogous 
to the scheme outlined in Ref. Here we have defined 
the single-particle matrix elements as 



i 



(B4) 
(B5) 



The matrix elements between Slater determinant and 
also those of single particle operators between different 
Slater determinants as used e.g. in Eq. (jA12[) are given 
by 



(2?|1^> (2?|2^> 



(2?|A^^> 



(B6) 



+ 



+ 



(l"l".,.|lg> 

(2?l".,.|lg) 



(l"|2g) 
(2?|2g) 



•• (KINS) 



+ 



{K\4) (l^|2g> 
(2^|lg> (2?|2g) 



(B7) 
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